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We study theoretically and numerically the conditions required for the appearance of a super¬ 
radiance transition in graphene. The electron properties of graphene are described in the single 
p^-orbital tight-binding approximation, in which the model is reduced to the effective two-level 
pseudo-spin 1/2 system. For each level we introduce the electron transfer rate of escape into a con¬ 
tinuum. We demonstrate that, under some conditions, the superradiance experiences the maximal 


in a wide THz frequency range. 

Qualitatively, the ST occurs when the reso¬ 
nances begin to overlap - the spacing between 
the resonances becomes of the order of the sum 
of half-widths of these resonances. With further 
overlapping of resonances, segregation of the 
eigenenergies takes place, depending on their 
decay widths. Namely, the wide superradiant 
eigenstates provide rapid and coherent decay of 
the initially populated state in the continuum. 
The subradient eigenstates, with narrow decay 
widths, survive for a relatively long time. Note, 
that not only the overlapping of resonances and 
the segregation of the eigenenergies are impor¬ 
tant for the occurrence of the superradiance, 
but the initial population of the system is also 
important. All these effects are described in 
details in Hflini- 

In this letter, we determine theoretically and 
numerically the conditions at which the ST oc¬ 
curs in a single-layer graphene material. We 
demonstrate how the ST in this system is re¬ 
lated to both the occurrence of an exceptional 
points (EP), when complex eigenvalues coin¬ 
cide, and to the overlapping of two resonances. 
We also show that, under some conditions on 
the parameters and initial populations, the 
maximal efficiency of the electron transfer (FT) 
into the sink is related to the ST. We compare 
the ST criterion based on the overlapping of 
resonances and on the occurrence of the EP. 
We 


quantum coherent escape to the continuum. 

PACS numbers: 42.50.Nn, 72.80.Vp, 03.65.Vf 
Keywords: Superradiance, graphene 


I. INTRODUCTION 


The superradiance transition (ST) was first 
described by Dicke in 1954 [l|. The ST is 
normally associated with a significant enhance¬ 
ment of the spontaneous radiation due to quan¬ 
tum coherent effects. Later it was demon¬ 
strated that the ST occurs in many quantum 
optical systems, nuclear systems (heavy nu¬ 
clei decay), nano- and bio-systems [J-Q. The 
ST usually occurs when the discrete (intrinsic) 
states of the system interact with the contin¬ 
uum spectra (sinks). Then, an adequate ap¬ 
proach for describing the eigenstates and the 
dynamics of the system can be based on an ef¬ 
fective non-Hermitian Hamiltonian for intrin¬ 
sic states In this case, the eigenener¬ 

gies of the non-Hermitian Hamiltonian become 
complex. Recently, a stady-state superradiant 
laser with less than one intracavity photon was 
demonstrated with rubidmm-87 atomic dipoles 
0. Large enhancement of Forster resonance 
energy transfer on graphene platforms was dis¬ 
cussed in 0 . In [13 , a superradiant plas- 
monic lasing with a giant gain at the palsmon 
modes in graphene was theoretically analyzed 
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FIG. 1: (Color online) Graphene two-dimensional 
lattice, with two atoms, A (red) and B (green), 
in a unit cell. Primitive unit vectors: a. 1^2 — 
(a/2)(v^, 3), b = (0, o), a « 1.42 A. 


II. DESCRIPTION OF THE MODEL 



FIG. 2: (Color online) The band structure of 
graphene. The valence and conduction bands touch 
each other in six DP points, where E(k.) = 0. 


The structure of a single atomic layer 
graphene can be described by the honeycomb 
lattice which consists of two triangular Bravais 
sublattices, represented in Fig. [T] by nonequiva¬ 
lent A (red) and B (g reen) carbon atoms, which 
create a unit cell [l2|, • Both sublattices have 

the periodic structures, and are shifted by a 
vector, b = (0,a), a ~ 1.42 A, which connects 
A and B atom in the unit cell. 

In a carbon atom, six electrons occupy the 
2s^, and 2p orbitals. From them, four 
valence electrons are responsible for structural 
and electronic properties. One of the valence 
electrons of each A and B atoms occupies the Pz 
orbital, which is orthogonal the graphene plane. 
The hybridization of these Pz orbitals pro¬ 
vides the formation of the 7r-bands in graphene. 
Then, the electron properties of graphene can 
be described within a single Pz orbital tight- 
binding (TB) model M- Using the TB approx¬ 
imation, one can show that the effective single¬ 
electron Hermitian Hamiltonian is reduced to 
the two-level pseudo-spin one-half system. The 
projections of pseudo-spin are associated with 
two sublattices. 

In graphene, the dispersion relation, ^(k) 
(where, k, is the wave vector), has some spe¬ 
cific properties: the Fermi level corresponds 
to ^(k) = 0, and the valence and conduction 
bands touch each other in the first Brillouin 
zone at six points. Each of these points pro¬ 
vides a “conical intersection” known also as the 
“diabolical point” (DP) [l3,[HI- (See Fig. O) 

According to 0 , in the k-representation the 
dynamics of the electron, in the vicinity of the 


DP, can be described by the following Hamilto¬ 
nian: 


Ho 


0 hvF{qx-iqy) 

hvF{qxEiqy) 0 


( 1 ) 


where, vf = 3\t\a/2h, is the Fermi veloc¬ 
ity. Here t ^ —2.8eV, is the hopping in¬ 
tegral for nearest neighbor atoms, A and B, 
with coordinates, and R^- We use the 
following notations: the state, |1) = (q): 
corresponds to the population of the sublat¬ 
tice A, and the state, |2) = ( ^ ), corresponds 
to the population of the sublattice B. The 
Hamiltonian, Ho^ has the eigenvalues, E± = 
EhvF\c[\- The corresponding eigenstates are: 
1+) = (l/v^)(e-*^/2|l) + e*^/2|2)) (conduction 
band), and |—) = (1/Ay^)(e“*'^/^|1) — e*'^/^|2)) 
(valence band), where p = arg(ga, + ig^). 

Suppose that each discrete state (related to 
the sublattice, A and B) is coupled to a con¬ 
tinuum - continuous energy band, which could 
originate due to impurities or other mecha¬ 
nisms. Suppose, that we allow an electron to 
tunnel to the continuum from the sublattices, 
A and B, with the ET rates, Fi and F 2 , cor¬ 
respondingly. Then, the quantum dynamics of 
the ET can be described by the following effec¬ 
tive non-Hermitian Hamiltonian, EL = EL — iW, 
where, EL^ is the dressed Hamiltonian, ELq, and 


We find. 





s — iV 
V 


u* 

—£ T iF 


( 2 ) 

( 3 ) 
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where, Aq = 5o - iFo, 5o = + ^ 2 , Fq = 

(Fi + F2)/2, V = 2hvF{qx + s = ei - 62 , 
and F = (Fi — F 2 )/ 2 . Here, 5^, is the renor¬ 
malized energy of the state, |n), which usually 
occurs in the non-Hermitian Hamiltonian ap¬ 
proach [H-lil. Note, that for pure graphene at 
the DP, the “electron masses”, 51,2 = 0, for 
both sublattices. However, due to the finite 
bandwidths, associated with the sinks, not only 
the ET rates, Fi ^25 appear, but also the small 
“effective electron masses”, 51 ^ 2 , may occur. 

Below, in numerical simulations, we choose 
h = 1. All energy-dimensional parameters 
are measured in ps“^ ~ 0.66 meV. Time is 
measured in ps. For example, if we choose, 
|q| = 10^ cm“^ and vp = 10^ cm/sec, we have: 
\V\ = 2 ps“^ « 1.32meV. 


III. SUPERRADIANCE TRANSITION 


The bi-orthogonal eigenstates of the effective 
non-Hermitian Hamiltonian, H, provide a con¬ 
venient basis in which the eigenvalue problem 
can be formulated and resolved. The solution 
of the eigenvalue problem, Hlu) = E\u) and 
{ulH = E{u\ (where, \u) and (fx|, are the right 
and left eigenvectors respectively), is given by. 


^ 1,2 = Ao/2 ± U/2, (4) 


where, U = \/|Hp + (e — iF)^. The right and 
the left eigenstates form the bi-orthonormal 
basis with the properties: (fii, 2 |'^ 2 ,i) = 0, 
(fil,2|'^l,2) = 1- 

Further, it is convenient to set, Ea = Sa — 
iTa = 1,2), where, £a = and = 

—is the value of the complex enemy (or 
the half-width of the resonance, a [3|-[6|). In 
Fig. m the values, Tq,, as the functions of Fi 
and F 2 , are presented. 

For Fi = F 2 , one has the standard en- 
er gy levels repulsion, A£ = Uq, where Uq = 
y/ 5 ^ + |Up, and A£ = £i — £2 denotes the en¬ 
ergy level spacing. When |F| ^ Uq, we obtain. 


A£ 

Ti 




T2 


£2 

2 


(5) 


Uq 


F. (6) 





(a) 





From here it follows that, while the decay FIG. 3: (Color online) Decay widths, (o = 1, 2), 

widths separate as F grows (a segregation effect as the functions of the rates, Fi and F 2 (|F| = 2): 

(H), the real energies, £i and £ 2 , are attracted ^ = 0; (c,d) £ = 2. 

to each other. 
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In the opposite limit, |r| ^ Qq, we find, 


AS ^ 

Ti« 


5 + 


2r2’ 


£1 

2 2r’ 


T 2 


£2 

2 


(7) 


2r' ^ ^ 


One can see that if Fi r 2 , the eigenstate, 
\ui), corresponding to ^ 1 , is the superradiant 
state, and the eigenstate \u 2 ), corresponding to 
^ 2 , is a long-living subradient state. In the 
opposite case, r 2 ^ Fi, the eigenstate, \ui), 
becomes long-living subradient state, and the 
eigenstate, \u 2 ), is a rapidly decaying superra¬ 
diant state. 

For two-level system, governed by the effec¬ 
tive non-Hermitian Hamiltonian, the quali¬ 
tative criterion of the overlapping of two reso¬ 
nances can be written as follows: Tq/AE ^ 1 , 
which, after some algebra, reduces to the equa¬ 
tion. 


r^ + r2r2-^^2r2-£2r2 = o. (9) 

When F 2 = 0 (or Fi = 0), the solution 
of this equation takes a simple form. Ft 9 = 

In Fig. m the half-widths of the resonances, 
Tq,, and the spacing, Af, vs Fi are depicted, 
for the choice of the FT rate, F 2 = 0. The ST 
occurs at the point, where the smallest half¬ 
width, T 2 , reaches its maximum value (Fig. 
m- One can see that the qualitative criterion 
of the overlapping of resonances is good enough 
for 5 > 2. Indeed, as Fig. |4] demonstrates, for 
6 = 2 (and \ V\ = 2 ), the ST occurs at the point 
Fi « 5.3, while the criterion of the overlapping 
of resonances yields, F^ 4.9. The error in the 
estimate of the ST is, « 7.5 %. 

In the opposite case, of the small values of 
5 ^ 1, as in graphene, the ST corresponds to 
the EP {s = 0), and the criterion of the over¬ 
lapping of resonances cannot be used. 


A. Superradiance transition in the 
vicinity of the exceptional point 

The effective non-Hermitian operator H has 
a singularity when the eigenvalues Ei and E 2 
coincide. To describe the behavior of the eigen¬ 
values in the vicinity of the degeneracy, it is 
convenient to consider the parameters: X = 
5RE, Y = and Z = e — iT. In the com¬ 
plex parameter space (X, T, Z), crossing points, 
known also as exceptional points (EP) [l 6 l-[^. 
are determined by equation, + = 0 . 



FIG. 4: (Color online) Left panel: (o = 1, 2) 

as the functions of the rate, Fi; Ti (blue) and T 2 
(red). Right panel: Dependence of the spacing, 
AS (blue), on Fi. Green line represents the sum, 
T = Ti + T 2 , vs Fi. Parameters: F 2 = 0, |F| = 
2, £ = 0 (solid curves); £ = 0.1, (dashed curves); 
£ = 0.5 (dotted curves); £ = 2 (dot-dashed curves); 
e = A (long-dashed curves). 




FIG. 5: (Color online) The behaviors of the real 
parts of the eigenvalues, Sa = 5R(Ea), in the vicin¬ 
ity of the EP, for F = 1. Left panel: £ = 0. The 
EP is represented by a circle, X^ -h = F^. Right 
panel: e: = 0.25. 


In Fig. O the real parts of the complex eigenen- 
ergies, S^^ are depicted. For 5 = 0 (left panel), 
the eigenvalues coalesce at the EP, represented 
by the circle, X^ + = F^. For 5 = 0.25 

(right panel), the gap between the valence and 
conduction bands occurs. 

Note, that in contrast to the case of the Her- 
mitian Hamiltonian, where the degeneracy is 
referred to a “conical intersection” at the DP 
and the coalescence of eigenvalues results in dif¬ 
ferent eigenvectors, at the EP the eigenvectors 
merge, forming a Jordan block. This leads to 
the violation of the normalization condition at 
the EP: {ui^ 2 \ui^ 2 ) = 0. 

Let us assume that the renormalization ener¬ 
gies are the same for both sublatticies, Si =52- 
Then, for |F| < |E|, we obtain for the spacing, 
AS = VX^ -1- — F^, and both eigenstates 

have the same widths. From here, it follows 
that the real parts of the eigenenergies. Si and 
^ 2 , attract each other as |F| grows. At the EP 
the levels cross, and the complex eigenenergies 
loose their analytical properties [j-Q . 

For |F| > |E|, the real parts of eigenenergies 
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become equal, and the decay widths of the res¬ 
onances separate, |Ti — T 2 I = — X‘^ — Y‘^ 

(see Fig. IH solid lines). In particular, for 
|r| > |y|, we obtain Ti = ri/2^0{{\V\/\T\)^) 
and T 2 =r 2/2 + 0((|y|/|r|)2). 

Thus, if Fi ^ r 2 , one has a superradiant 
eigenstate with F^i, and the subradient long- 
living eigenstate with E 2 . In the opposite case, 
r 2 ^ Fi, the eigenstate with Ei becomes sub¬ 
radient long-living, and the eigenstate with E 2 
becomes a superradiant rapidly decaying eigen¬ 
state. 


IV. THE ELECTRON TRANSFER 
INTO SINKS 

The FT quantum dynamics can be described 
by the Liouville-von Neumann equation, 

ip = i[p,n]- {W,p}, ( 10 ) 

where {W, p} = Wp + pW. We define the ET 
efficiency of tunneling into n-th sink as min, 

r?„(i)=r„/ pnn{T)dT, n = 1,2. (11) 

Jo 


So, the efficiency, pi, decreases for large values 
of Fi. 

As Fig. [6] shows, even if only the conduction 
band (dot-dashed curves) or the valence band 
(dotted curves) of the Hamiltonian ([T]), is ini¬ 
tially populated, the maximum in pi can still 
be observed. Evidently, when the site state, |1), 
is initially populated, the maximum in pi does 
not occur. 

In Eig. [71 the efficiencies, pi and p 2 , are 
demonstrated as functions of the ET rates, Fi 
and r 2 , for initially population of the sublat¬ 
tice B, and for relatively small e = 0.1ps“^ ~ 
0.066meV. 



r 

1 


Employing Eq. m, one can show that 
the following normalization condition is sat¬ 
isfied, Trp(t) -I- p(t) = 1. Here, p(t) = 

/q Tr{>V,p(r)}(ir = pi(t) +P2(0, is the total 
ET efficiency of tunneling into both sinks. 

In Fig. ini the efficiencies, pi (blue curves) 
and r ]2 (red curves), are shown as functions of 
Fi, for fixed values of r 2 = 1 and e = 0.1. 
The time of evolution was chosen, r = 10 ps. 
The choice of the time-interval, r, was based on 
the numerical experiments which demonstrated 
that during this time the efficiencies of sinks 
reached their asymptotic values. 

If the sublattice B is initially populated, the 
efficiency pi (blue solid line) experiences a max¬ 
imum at Fi ^ 2. This maximum is associated 
with the superradiance into first sink connected 
with the sublattice A, and it can be qualita¬ 
tively explained as follows. The initial popula¬ 
tion of the sublattice B corresponds to the site 
state, 12). When Fi ^ r 2 , the ET mainly oc¬ 
curs into the second sink. So, p 2 is close to 1 
and Pi ^ 1 in this region. When Fi increases, 
the ET occurs with maximum in rj 2 (blue curve) 
and minimum in pi (red curve). When Fi 
further increases, the superradiant eigenstate 
(with large decay width, Fi) is mainly local¬ 
ized on the site state, |1), which only weakly 
overlaps with the initially populated state |2). 


FIG. 6: (Color online) The efficiencies of sinks, pi 
(blue) and p 2 (red) vs Fi, in the absence of noise. 
Initial conditions: pii(O) = 0, pi2(0) = P2i(0) = 0 
P22(0) = 1 (solid curves); pii(O) = 1, pi2(0) = 
P2i(0) = 0 p22(0) = 0 (dashed curves); pii(O) = 
0.5, pi2(0) = P2i(0) = 0.5 p22(0) = 0.5 (dot-dashed 
curves); pii(O) = 0.5, pi2(0) = p2i(0) = -0.5 
P22(0) = 0.5 (dotted curves). Parameters: r 2 = 1, 
A = 2, y = 0, £ = 0.1. Time of evolution: 
T = 10 ps. 



FIG. 7: (Color online) The efficiencies of the sinks 
vs Fi and F 2 , in the absence of noise. Initial con¬ 
ditions: Pll(o) = Pl2(0) = P2l(0) = 0 , P22(0) = 1. 
Parameters: |y| = 2, e = 0.1. Time of evolution is: 
T = 10 ps. 


Using the results obtained in [^, we derive 
the analytical expressions for the efficiencies, pi 
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and ? 72 - For illustrative purposes, in what fol¬ 
lows, we restrict ourselves by some simple cases. 

Let us assume that initially the sublattice B is 
populated. Then, the efficiency of trapping the 
electron into the first sink is given by, T] 

Vi (t) = Vo — (^(Fo cosh + ^2 sinh f^ 2 ^) 

—C(rocosf^it — fti sinf^it)), (12) 


01 23456789 10 

r 

1 


FIG. 9: (Color online) The efficiencies of sinks with 
noise, 771 (blue) and 772 (red) vs Fi: d = 0 (solid 
curves), d = 5 (dashed curves), d = 10 (dotted 
lines), d = 20 (dot-dashed curves). Parameters: 
7 = 10 , F 2 = 1 , |y| = 2 , d = di - d 2 , £ = 0 . 
Initial conditions: pii(O) = 0, pi2(0) = P2i(0) = 0 
^ 22 ( 0 ) = 1. Time of evolution: r = 10 ps. 


where, Cl 2 = and 


\V\^ToT, 

2iTl + nj)(rl-niy 

(13) 

p Vpri 

2(e2 + el)(r2-ei)’ 

(14) 

2(^1 + ^2)(Fo + ^ 1 ) 

(15) 




the maximum of the efficiency, 770, in the first 
sink depending on both parameters. The most 
pronounced maximum of 770 corresponds to the 
graphene case - absence of the effective electron 
mass, 5 = 0. 


V. NOISE-ASSISTED ELECTRON 
TRANSFER 


In the presence of noisy environment, the 
evolution of the system can be described by the 
following effective non-Hermitian Hamiltonian 


FIG. 8 : (Color online) The asymptotic efficiency of 
the 1 st sink, 771 ( 00 ) = 770 , vs Fi and 5 , in the ab¬ 
sence of noise. Initial conditions: pii(O) = pi 2 ( 0 ) = 
P2i(0) = 0, p22(0) = 1. Parameters: |y| = 2, 
F 2 = L 

For I Ft I ^ 1, the asymptotic behaviors of 
771^2 are given by, 

7ii{t) ^rjo- (16) 

ri2{t) ^ 1 - rio + (17) 


ntot = n-iW + V{t), (19) 

where, V{t) = J2m,n^rnn{t)\m){n\, m,n = 
1,2, and Xmn{t) describes the influence of the 
noisy environment. In what follows, we re¬ 
strict ourselves with the diagonal noise, writing 
Xmn {t) = Xn^mn^{t)^ where is the coupling 
constant at site, n, and ^(t) is the random tele- 
graph process (RTP) with the properties Hi- 


where, 

4{nj + nl){ro-n2)' ^ ^ 

For fixed F 2 , the behavior of the asymptotic 
value, 770, as the function of Fi and 5 is demon¬ 
strated in Fig. [HI One can see the details of 


where, x(t—t') = I, is the correlation 

function of noise. 

The evolution of the average components of 
the density matrix is described by the following 
closed system of ordinary differential equations 
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ilii: 

j^{p)=i[{p),H]-{W,{p)}-iB{p^), ( 21 ) 

|(p«) =i[(/),H]-{W,(p«) 

-iB{p)-2j{p^), (22) 

where (/?«) = {^p)/a, B = 'Em,ni^rn - 

dn)\m){n\^ and we set dm = A^cr. The aver¬ 
age, (...), is taken over the RTF. 

In Fig. [9l the influence of noise on the ET 
rate is demonstrated, for the effective electron 
mass, £ = 0 , and for fixed parameters: r 2 = 1 , 
7 = 10 and \V\ = 2 . The amplitude of noise 
is characterized by the parameter, d = di — d 2 . 
As one can see from Fig. [9l the ET is stable 
in the sense that when d 7 ^ 0 , the maximum 
in r]i still exist, but is less pronounced, when d 
increases. 


VI. CONCLUSION 

We studied theoretically and numerically a 
possibility to observe a superradiance transi¬ 
tion in a single-layer graphene material. The 
analysis was performed in the tight-binding ap¬ 
proximation and at the vicinity of the zero 
gap between the conduction and the valence 
bands. For each triangle sublattice we intro¬ 
duced the electron transfer rates in their con¬ 
tinuum (sinks). So, the total model, in the k- 
representation, is reduced to the two-level sys¬ 
tem governed by the non-Hermitian Hamilto¬ 
nian. 

In the graphene case, the superradiance tran¬ 
sition occurs in the vicinity of the exceptional 
point, where two complex eigenenergies coin¬ 
cide. Relatively far from the exceptional point. 


a qualitative interpretation of the superradi¬ 
ance transition is associated with the overlap¬ 
ping of two resonances: when the sum of half¬ 
widths (decay rates) of resonances is of the or¬ 
der of the spacing between them (distance be¬ 
tween the real parts of the complex eigenener¬ 
gies). When one of the electron transfer rates 
further increases, a segregation of the reso¬ 
nances occurs. Namely, one decay rate, corre¬ 
sponding to the superradiant eigenstate, grows 
and the second one, corresponding to the sub- 
radient eigenstate, becomes limited from above. 

The superradiant state is usually responsible 
for rapid coherent electron transfer into a con¬ 
tinuum (sink). Under some conditions on pa¬ 
rameters and initial population, which are ana¬ 
lyzed in the paper, the efficiency of the ET into 
a sink has a maximum at the particular value 
of the electron transfer rate. 

In order to observe the described superradi¬ 
ance in the graphene material, some additional 
experimental issues must be resolved. In par¬ 
ticular, it is required to realize and control two 
electron transfer rates for two sublattices. It is 
also required to prepare the system in the par¬ 
ticular initial states, say, to populate initially 
one of the sublattices. The experimental obser¬ 
vation of the discussed superradiance transition 
will be important for both, better understand¬ 
ing of the fundamental properties and for con¬ 
trolling of the graphene based material. 
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